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Abstract. Cornerstone molecules (CO, H2CO, CH3OH, HCN, HNC, CN, CS, SO) were observed toward seven sub-millimeter 
bright sources in the Orion molecular cloud in order to quantify the range of conditions for which individual molecular line 
^ ' tracers provide physical and chemical information. Five of the sources observed were protostellar, ranging in energetics from 

^ : 1 - 500 Lq, while the other two sources were located at a shock front and within a photodissociation region (FDR). 

Statistical equilibrium calculations were used to deduce from the measured line strengths the physical conditions within each 
source and the abundance of each molecule. In all cases except the shock and the FDR, the abundance of CO with respect 
to H2 appears significantly below (factor of ten) the general molecular cloud value of 10"*. Formaldehyde measurements 
were used to estimate a mean temperature and density for the gas in each source. Evidence was found for trends between 
I the derived abundance of CO, H2CO, CH3OH, and CS and the energetics of the source, with hotter sources having higher 

abundances. Determining whether this is due to a linear progression of abundance with temperature or sharp jumps at particular 
temperatures will require more detailed modeling. The observed methanol transitions require high temperatures (T > 50 K), 
and thus energetic sources, within all but one of the observed protostellar sources. The same conclusion is obtained from 
^ ■ observations of the CS 7-6 transition. Analysis of the HCN and HNC 4-3 transitions provides further support for high densities 

i-^ ' n > 10' cm"^ in all the protostellar sources. 

, The shape of the CO 3-2 line profile provides evidence for internal energetic events (outflows) in all but one of the protostellar 

sources, and shows an extreme kinematic signature in the shock region. In general, the CO line and its isotopes do not signif- 
icantly contaminate the SSOfiin broadband flux (less than 10%); however, in the shock region the CO lines alone account for 
more than two thirds of the measured sub-millimeter flux. In the energetic sources, the combined flux from all other measured 
' molecular lines provides up to an additional few percent of line contamination. 



Key words, astrochemistry ~ ISM:molecules(CO,H2CO,CH30H,HCN) - ISM:individual(Orion A:OMCl,OMC2,OMC3 - 
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1. INTRODUCTION 

The study of structure within the star-forming regions of 
molecular clouds has benefitted significantly from observations 
of many molecular lines, each tracing specific chemical and 
physical conditions. Most studies either focus on surveys using 
particular tracers, such as carbon monoxide, N2H^, or ammo- 
nia, in order to map the column density and kinematics of the 
gas (Myers 1999), or have focused on individual sources, pro- 
ducing detailed molecular catalogues in several to dozens of 
species designed to constrain the physical and chemical mor- 
phology of the region (van Dishoeck 8c Blake 1998). 
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Several molecules have emerged as particularly good di- 
agnostics of the conditions and chemistry near young stellar 
objects. For example, H2CO has many lines that are readily 
observed at sub-millimeter wavelengths, and whose ratios are 
either good temperature or density tracers (e.g., Hurt, Barsony, 
& Wooten 1996, Mangum, Wooten, & Barsony 1999, Mitchell 
et al. 2001). Analysis of dust continuum and line emission has 
shown that temperature and density gradients exist across the 
protostellar envelopes, with temperatures varying from the in- 
ner to the outer region from > 100 K to less than 20 K and den- 
sities from > 10^ cm"^ to ~ lO'* cm"^ (e.g., van der Tak et al. 
2000a, Shirley, Evans, & Rawlings 2002, J0rgensen, Schoier, & 
van Dishoeck 2002). The chemistry responds to these changes. 
Molecules freeze-out onto the grains in the cold outer parts of 
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the clouds, where they can form new species through grain- 
surface reactions. The composition of icy mantles can be deter- 
mined through infrared absorption, with species such as H2O, 
CO, CO2 and CH3OH known to have high ice abundances (e.g., 
Gibb et al. 2000). In the inner region close to the protostar, the 
grains are heated and the ices are observed to evaporate back 
into the gas (e.g., van Dishoeck & Helmich 1996; Boonman & 
van Dishoeck 2003). High temperature gas-phase reactions be- 
tween evaporated species can subsequently lead to high abun- 
dances of complex organic molecules observed in high-mass 
hot cores (e.g., Rodgers & Charnley 2001). Depending on the 
evolutionary state of the source, different chemical characteris- 
tics become more prominent. 

In this paper we consider a selection of cornerstone 
molecules and study a variety of independent locations within 
the Integral Shaped Filament (ISF) in Orion A, along with NGC 
207 IIR in Orion B, in order to quantify the range of condi- 
tions for which individual molecular line tracers provide phys- 
ical or chemical information. This is primarily a morphologi- 
cal study, comparing the differences in molecular line emission 
across sources of varying physical conditions using simplified 
equilibrium modeling, in order to search for obvious diagnos- 
tic features within the data. The list of sources observed in 
Orion A was compiled from a sensitive dust continuum survey 
of the ISF (Johnstone & Bally 1999) and includes both highly 
enshrouded sub-millimeter bright sources, possibly protostel- 
lar, through more evolved protostars, a bright PDR knot, and 
a shock front. Many of the young stellar objects are of inter- 
mediate mass and thus the observations bracket sub-millimeter 
studies of high-mass (e.g. Hatchell et al. 1998; van der Tak 
et al. 2000a) and low-mass (e.g. Blake et al. 1994, 1995; van 
Dishoeck et al. 1995; Buckle & Fuller 2002; J0rgensen et 
al. 2002) star-forming regions. The selection of cornerstone 
molecules and molecular transitions was chosen based on pre- 
vious experience in detailed studies of individual objects. 

2. The Orion Integral Shaped Filament 

The ISF (Bally et al. 1987) in the northern portion of the Orion 
A molecular cloud is the nearest site with active high-mass 
through low-mass star formation. The region is associated with 
the Orion Nebula and the Trapezium Cluster of 700 young 
stars (Hillenbrand 1997; Hillenbrand & Hartmann 1998), con- 
tains the OMCl cloud core immediately behind the Nebula, 
and two other extensively studied star-forming cores, OMC2 
and OMC3 located about 15' and 25' to the north (Castets 
& Langer 1994). To the south the region becomes much qui- 
eter with OMC4 (Johnstone & Bally 1999) showing no bright 
sub-millimeter sources. The ISF has spawned several thousand 
young stars in the past few million years and contains dozens 
(possibly hundreds) of embedded young stellar objects (Chini 
et al. 1997, Johnstone & Bally 1999) which power dozens of 
molecular outflows, Herbig-Haro objects (Reipurth, Bally, & 
Devine 1997), and molecular hydrogen emitting shocks and 
jets (Yu, Bally, (& Devine 1997). The abundance of star for- 
mation, at various stages of evolution, within a single, nearby, 
molecular cloud provides a unique setting for conducting tests 
of astrochemical diagnostics. 



2.1. The Individual Sources 

The sources chosen for this study were primarily drawn from 
the 450 and 850 /urn dust emission maps of the ISF (Johnstone 
& Bally 1999; Fig.Q). Supplemental information on the prop- 
erties of the sources comes from measurements at 350/im 
(Lis et al. 1998), 1300 //m (Chini et al. 1997), 3.6 cm VLA 
data (Reipurth, Rodriguez, & Chini 1999), and the IRAS point 
source catalogue. Additional 2 yum molecular hydrogen obser- 
vations (Yu et al. 1997) of the ISF provide evidence of local 
heating expected within shocks, and delineate the location of 
molecular jets. The seven sources were chosen to cover a wide 
selection of physical conditions in order to test the hypothesis 
that molecular line observations provide sensitive diagnostics. 

The coldest sub-millimeter source in our sample is MMS6, 
T^i ~ 15K as measured by the sub-millimeter fluxes from 
350 to 1300 yum (Lis et al. 2001). MMS6 is also the bright- 
est 850 jum source within the OMC3 region of the ISF, with a 
peak flux of 7.4 Jy within a 15 arcsecond beam (Johnstone & 
Bally 1999). The total mass within the envelope is estimated 
to be 36 Mq (Chini et al. 1997). There is no corresponding 
IRAS point source, although the background brightness of the 
region makes point source detection difficult, suggesting that 
there is no strong internal heating source. The sub-millimeter 
luminosity, energy emitted at wavelengths greater than 350 fj.m, 
is 1.2 Lq; however, the lack of far-infrared emission measure- 
ments provides only an upper limit on the total bolometric lu- 
minosity of < 60 Lq (Chini et al. 1997). Despite the very cold 
appearance of MMS6 in the sub-millimeter, and thus the possi- 
bility that it is a pre-stellar clump heated only from the outside 
by the strong interstellar radiation field (ISRF) in Orion, MMS6 
has been classified as a Class protostar (Chini et al. 1997). Yu 
et al. (1997) observed knots of H2 directly north and south of 
MMS6 which are interpreted as bow shocks within a molecular 
outflow from the continuum source. Chini et al. (1997) report 
substantial high-velocity red-shifted emission in '^CO centered 
on MMS6, although no corresponding blue-shifted lobe was 
observed. More recent observations by Aso et al. (2000) as- 
sociate the CO outflow with a different sub-millimeter source. 
VLA observations at 3.6 cm (Reipurth et al. 1999) locate a ra- 
dio source coincident with MMS6. Such coincidences are com- 
mon among Class sources and have been interpreted as ther- 
mal Bremsstrahlung radiation produced in shocks within the 
inner part of the molecular outflow. Thus, while MMS6 is pos- 
sibly a pre-collapse clump, it is likely that the sub-millimeter 
source enshrouds a very young Class protostar 

Despite being three times fainter than MMS6 at 850yum, 
MMS9 appears to be the source of the most powerful out- 
flow in the OMC2/3 region (Yu et al. 2000). The spectral in- 
dex for MMS9 is less steep than observed in MMS6, suggest- 
ing a warmer temperature in the dust envelope, T^i ~ 20 K 
(Chini et al. 1997). The total envelope mass is estimated to be 
10 Mq, considerably lower than the envelope mass of MMS6. 
The IRAS flux limits constrain the luminosity of the embedded 
source to be < 94 Lq (Chini et al. 1997). Despite having no 
coincident IRAS point source, MMS9 is expected to contain 
a well-advanced Class protostar. Reipurth et al. (1999) ob- 
served 3.6 cm continuum emission at this location. The rising 
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infrared background toward MMS9 may be responsible for the 
lack of an IRAS point source detection. 

The only IRAS point source within the OMC3 region is 
IRAS 05329-0505, coincident with MMS7. From the wave- 
length of the peak emission, /li^^x ~ 60;um, the warm com- 
ponent of the dust can be estimated, ~ 50 K. In the sub- 
millimeter MMS7 appears quite similar to MMS9, although the 
envelope temperature is somewhat higher ~ 26 K (Chini et 
al. 1997). The total envelope mass is ~ 8M0, comparable to 
the envelope around MMS9, and the bolometric luminosity is 
~ 76 Lq. Reipurth et al. (1999) observed 3.6 cm continuum 
emission at this location, providing further evidence for an em- 
bedded Class source. MMS7 appears to be in a somewhat 
more advanced stage of star formation, powering a giant bipo- 
lar Herbig-Haro flow (Reipurth et al. 1997) and harboring both 
a central FU Ori-like source (Reipurth & Aspin 1997) and an 
optical reflection nebula (Wolstencroft et al. 1986). 

Further south along the ridge, near the center of the OMC2 
molecular core, lies the far infrared source FIR4 (Mezger, 
Zylka, & Wink 1990), also observed as an IRAS point source. 
The region around FIR4 harbors a cluster of near-infrared 
sources (Johnson et al. 1990) although none are explicitly co- 
incident with the far-infrared source. FIR4 is the brightest sub- 
millimeter source in the 0MC2 region and contains a VLA ra- 
dio source (Reipurth et al. 1999), leading to identification as 
a Class object. Its proximity to the near IR cluster suggests, 
however, that it might be a more evolved envelope. The pres- 
ence of far-infrared emission provides evidence for warm dust 
Td ~ 40 K (Mezger et al. 1990), although the sub-milHmeter 
emission traces the much colder outer envelope ~ 15 K. The 
total envelope mass is ~ 34 Mq and the luminosity is ~ 400 Lq 
(Mezger a/. 1990). 

The final protostellar source in our sample is NGC2071IR 
in the Orion B North molecular cloud (Johnstone et al. 2001; 
Motte et al. 2001; Harvey et al. 1979). Significant far-infrared 
emission from the source (Butner et al. 1990) provides evi- 
dence of an internal source and a warm dust > 50 K com- 
ponent providing a reprocessed luminosity ~ 500 Lq. The total 
envelope mass, derived from the sub-millimeter dust emission, 
is 12 - 30 Mo (Johnstone et al. 2001). NGC2071 contains one 
of the most powerful outflows known (Bally & Lada 1983). 

Along with the five protostellar candidates, two sub- 
millimeter bright regions which are not associated with indi- 
vidual protostars were also observed. A sub-millimeter bright 
location, PDRl, along the warm photodissociation region (the 
Orion Bar) separating the HII region M42 (the Orion Nebula) 
from the Orion A molecular cloud was included to monitor 
changes in the molecular line diagnostics in regions of en- 
hanced UV fluxes and higher excitation temperatures (see also 
Hogerheijde, Jansen, & van Dishoeck 1995). Also, a strong 
shock front, SKI, located via bright 2/im H2 emission (Yu 
et al. 1997), was included to test for shock diagnostics. The 
shocked region is also visible in the spectral index map of 
the ISF (Johnstone & Bally 1999) as a region with anoma- 
lously shallow slope, requiring either a lower temperature than 
the surrounding material, a shallower dust emissivity profile, 
or contamination of the 850//m continuum measurement by 
molecular fines. 



A summary of the physical properties of the seven sources 
is presented in Tabled The sources are ordered from the cold- 
est, least energetic, to the warmest, most energetic, with the 
PDR and shock front noted separately. 

3. Observations 

3.1. Choice of Molecules 

Eight different molecules and their isotopes have been ob- 
served in order to find useful diagnostics to study physical 
and chemical differences between the sources in our sample. 
Among these are carbon-, oxygen-, nitrogen-, and sulphur- 
bearing species, allowing for a study of differentiation in the 
chemical abundances. Abundant CO is observed toward all 
sources. Nevertheless, it is a useful molecule to study deple- 
tion effects, which are expected to be largest for the coldest, 
less evolved sources. In those sources where CO is not too 
depleted, it quickly becomes optically thick. In that case, the 
optically thin isotope C'^O can be used to derive H2 column 
densities. 

Formaldehyde, H2CO, is a slightly asymmetric rotor 
molecule, providing it with many of the temperature and den- 
sity diagnostic features of purely symmetric rotors. Mangum 
& Wootten (1993) and van Dishoeck et al. (1993) fist many 
of the strengths of this molecule as a probe of the physical 
conditions in dense molecular clouds. Methanol, CH3OH, is a 
complex molecule with many transitions in the sub-millimeter. 
Since these fines are concentrated in complex bands spanning 
relatively small wavelength ranges, this allows for efficient ob- 
servations even at relatively high spectral resolution. Methanol 
is also chemically important, since it has been detected both in 
the gas-phase and on grains (e.g. Dartois et al. 1999). The high 
gas-phase abundances of CH3OH most likely result from evap- 
oration off of grains (Charnley, Tielens, & Millar 1992; van der 
Tak, van Dishoeck, & CaselH 2000b; Schoier et al. 2002). 

One of the most abundant nitrogen-bearing species is HCN. 
Chemical models predict that the HCN abundances increase 
with temperature (Rodgers & Charnley 2001), thus provid- 
ing a useful evolutionary indicator. Observations of massive 
star-forming regions confirm this (Lahuis & van Dishoeck 
2000), although the situation for low-mass sources is less clear 
(J0rgensen et al. in prep.). Schilke et al. (1992) show that the 
abundance ratio of HCN and HNC depends on the density and 
kinetic temperature of the gas, providing a physical diagnostic. 
The chemical equilibrium of HCN, HNC, and CN also depends 
on the hardness of the radiation field. 

Finally, sulphur-bearing species have been associated with 
shocks and/or outflows (e.g. Orion; Blake et al. 1987). 
However, recent studies of sulphur-bearing species in massive 
star-forming regions other than Orion cannot distinguish be- 
tween an origin in the circumstellar envelope or the outflow 
(Keane et al. 2001a; van der Tak et al. 2003). Thus, it is still 
debatable whether sulphur-bearing species are good diagnos- 
tics of outflows. On the other hand, CS transitions have been 
used to derive detailed density structure within the molecular 
envelopes of star-forming regions, making this species a useful 
indicator of physical properties (van der Tak et al. 2000a). The 
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H2S/SO2 ratio has also been proposed as a chemical clock of 
hot core regions (Charnley 1997), although observational re- 
sults are still ambiguous (Hatchell et al. 1998; van der Tak et 
al. 2003). 

3.2. Sub-millimeter Observations 

The molecular line observations were taken at the 15-meter 
James Clerk Maxwell telescope (JCMT)' on Mauna Kea in 
Hawaii. The data were obtained during three separate observ- 
ing runs: January 1999, September 1999, and December 2000. 
Table 121 provides details of the instruments and settings used 
during each run, along with typical receiver conditions. For 
each source, Table|3]lists the molecular lines observed and the 
results of fitting one or two Gaussian components to the inte- 
grated line intensity, the full width at half maximum, and the 
peak intensity. 

The data were obtained in a beam switching mode, with 
a 10 arcminute switch in right ascension. The integrated line 
data were transfered from the antenna temperature scale to 
the main-beam brightness scale r^b by dividing by the main- 
beam efficiency ?/mb- For the A band receiver (210-270 GHz) 
the measured efficiency is ?7mb = 0.69, while for the B band 
receiver (330-370 GHz) rj^b = 0.60. The relevant beam sizes 
for the JCMT are 21 and 14 arcseconds at respectively 230 and 
345 GHz. The velocity resolution ranged from 0.2 - 0.5 km s ' . 
The rms noise after a typical integration time of twenty minutes 
was 0.30 K per channel for A band observations and 0. 15 K per 
channel for B band observations. Further details on the JCMT 
and the heterodyne receivers can be found on the JCMT home- 
page^. 

3.3. Representative Line Profiles 

Many important differences between the sources are readily ob- 
served simply by looking at the spectra. The extremely strong, 
and almost always self-absorbed, CO 3-2 transition reveals in- 
formation on the gas kinematics within each source (Fig. |2j. 
Single component Gaussian fits to the profiles are only pos- 
sible for MMS6 (line width 4kms-') and PDRl (3kms-'). 
MMS9 requires both a narrow component (6 km s ' ), providing 
two thirds of the integrated intensity, and a broad component 
(15 km s"'). FIR4 appears to be a scaled-up version of MMS9, 
with a narrow component (9kms ') containing three quarters 
of the integrated intensity and a broad component (21 kms '). 
Both the shock region and NGC207 1 are poorly fit by Gaussian 
components and are dominated by the extended wings. 

Eight of the methanol lines are observed in a single spectral 
setting and can be plotted for each source (see examples in Fig. 
|3}. These spectra also contain the sole observed SO transition. 
In agreement with the results found by van der Tak et al. (2003) 
there is an apparent lack of correlation between the strength 

' The JCMT is operated by the Joint Astronomy Centre on behalf 
of the Particle Physics and Astronomy Research Council of the UK, 
the Netherlands Organization for Scientific Research, and the National 
Research Council of Canada. 

^ _http : //www ■ j ach . Hawaii . edu/ J ACpublic/JCMT| 



of the SO line and the methanol lines. As well, differences in 
the relative strengths of the methanol lines are observable be- 
tween MMS6 and FIR4, indicating the more extreme temper- 
ature and density conditions in the energetic infrared source. 
Finally, the FIR4 methanol spectra show evidence for two com- 
ponents, narrow and broad, similar to the CO spectra. 

The HCN and HNC 4-3 profiles are also observed in a sin- 
gle spectral setting (Fig.|^. A range of relative line strengths is 
seen in the spectra; however, optical depth considerations make 
simple abundance determinations more complicated (see Sect. 
14.5. 1> . The substantial optical depth of the transitions allows 
for more detailed determinations of the source kinematics. In 
all sources except PDRl, the HCN line profile is broader than 
the HNC line profile or has evidence for extended line wings. 

4. Molecular Line Survey 

4.1. Statistical Equilibrium Calculations 

The analysis of the data is done with statistical equilibrium 
calculations, using the numerical code RADEX^. The adopted 
method is described in detail by Jansen, van Dishoeck, & Blake 
(1994), van Dishoeck et al. (1993), and Jansen et al. (1995). 
This method assumes that the source is a homogeneous spheri- 
cal region. The decoupling of radiative transfer and level popu- 
lations is done with the escape probability method. The optical 
depth T at the line center is proportional to the ratio of the total 
column density and the line width of the molecule. The optical 
depth and the line intensity depend only on this ratio. All vari- 
ables can be calculated, using a geometrical assumption that 
they are constant throughout the region. The radiative excita- 
tion of rotational lines in the sub -millimeter region is usually 
dominated by the 2.73 K cosmic background radiation, which 
peaks at 1.871 mm. Therefore, this external radiation field is in- 
cluded in the calculations. No central radiation field is included 
nor infrared pumping by dust. Since the sources have rather low 
dust temperatures, this latter effect is expected to be negligible. 

The kinetic temperature r^n and the density «(H2) for the 
homogeneous region can be estimated from line ratio calcu- 
lations where available. The column density of molecule X, 
N{X), can then be estimated using the temperature and the den- 
sity as found from line ratio calculations, such that the line 
width, integrated line strength, and main beam temperature of 
the emission lines can be fitted. 

Due to the simplistic treatment of the envelope structure 
the derived abundances are only roughly determined. Rising 
temperature and density gradients toward the envelope centers 
are expected due to heating by the central protostellar sources 
and these conditions can influence the observed line strengths 
significantly. Additionally, these Orion sources reside within 
an enhanced external radiation field, due to the proximity of 
many nearby massive stars (Li, Goldsmith, & Menton 2003). 
Thus, the outer envelopes are kept warm, ~ 30 K (Johnstone 
& Bally 1999). ModeHng of such envelopes (Li et al. 2003, 
J0rgensen et al. in prep.) reveal an inward decreasing tempera- 
ture gradient near the surface, encompassing much of the mass 

^ This program was originally written by J.H. Black. Here a more 
recent version by D.J. Jansen is used. 
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of the envelope, coupled to an inward rising temperature gradi- 
ent once the central heating source dominates. 

The large, externally warmed, envelopes surrounding the 
protostellar sources in Orion should allow for reasonable abun- 
dance estimates from the observed line intensities using statis- 
tical equilibrium calculations, except for those transitions re- 
quiring either very high densities n > 10^ cm"^ or high temper- 
atures r > 50 K. As a guide to the accuracy of the calculations 
presented in this paper, detailed modeling of the protostellar 
source IRAS 16293-2422, including careful derivation of the 
physical structure, has shown that abundance estimates from 
this simpler procedure are typically only accurate to within a 
factor of two or three (Schoier et al. 2002). 

4.2. Isotopes of Carbon Monoxide 

Within star-forming regions of molecular clouds, comparison 
between the line strengths of low-lying CO transitions and the 
intensity of the sub-millimeter dust continuum often suggest 
that there is a depletion of CO with respect to the dust, and 
by extension the molecular hydrogen abundance (see for ex- 
ample Mitchell et al. 2001, J0rgensen et al. 2002, Bergin et al. 
2002, Bacmann et al. 2002). The most common explanation for 
this observational discrepancy between CO line strength and 
dust continuum emission is that the CO has depleted through 
freeze-out onto the grains. The time scale for depletion scales 
inversely with the density of the region, and is shorter than 10^ 
yr for densities above a few xlO^ cm"^. Recent laboratory ex- 
periments on CO and H2O ices show that CO will remain com- 
pletely frozen out if the grain temperatures stay below ~20 K. 
Above 20 K, CO starts to evaporate back into the gas phase in a 
comphcated pattern depending on the ice structure, with a frac- 
tion of the CO remaining trapped in the H2O ice matrix until 
the entire mantle evaporates at ~100 K (Collings et al. 2003, 
Fraser et al. 2001). 

4.2.1 . Line strengths 

Table 0] details the observed line strengths of CO, '^CO, and 
C'^O 3-2 transitions obtained at the JCMT, along with an es- 
timate of the column of carbon monoxide along the line-of- 
sight. Using RADEX, the column is deduced from the inte- 
grated line strength of the C'^O 3-2 line, assuming that the 
'^CO peak temperature provides a reasonable estimate of the 
kinetic temperature of the gas. The RADEX calculations for 
CO include energy levels up to 7 = 25 and use collisional rate 
coefficients by Schinke et al. (1985). The results are not signif- 
icantly changed if the new rate coefficients by Flower (2001) 
are adopted (see J0rgensen et al. 2002). For NGC2071, where 
the '^CO line was not observed, the temperature was assumed 
to be 30 K, similar to the temperature of FIR4. 

The cosmic abundance ratio for isotopes of Oxygen, 
['^0]/['^0]~ 1800 (Wilson & Rood 1994), was used to convert 
the calculated C'^O column to an equivalent column of CO. 
For the expected conditions within the sub-millimeter clumps, 
n > 10"* cm"^, there is only a slight density dependence on the 
calculated column of the gas. In most cases the CO peak tem- 



perature was larger than the '^CO temperature; however, the 
optical thickness of CO is extreme and hence the line profile 
may be dominated by a small fraction of hot gas associated with 
each source. There is still uncertainty comparing the thicker 
'^CO with the thinner C'^O. We estimate the reliability of the 
computed physical parameters by comparing the RADEX com- 
puted optical depth of the '^CO line, Tmod, with the observation- 
ally determined optical depth Tobs, measured by comparing the 
peak temperature of the C'^O against the '^CO peak tempera- 
ture and assuming [i3CO]/[Ci^O]= 25 (Wilson & Rood 1994). 
This calculation confirms that the '^CO line is optically thick 
in all cases (tuco > 1-3). While there are deviations between 
the two determinations of the optical depth, the order unity dis- 
persion provides an equivalent estimate of the reliability of the 
measured column density of CO. 

4.2.2. Determination of CO abundance 

The column density of molecular hydrogen within each en- 
velope is best determined from the sub-millimeter flux mea- 
surements. Assuming that the 850;um emission observed by 
Johnstone & Bally (1999) and Johnstone et al. (2001) traces 
the bulk of the envelope and that the dust emission is optically 
thin, the mean column density through the center of each clump 
is related to the peak brightness: 



B850(Td)Af850A'mH, ' 

where Qbm is the effective solid angle of the SCUBA beam (in 
steradians), B^=,q(T^) is the Planck function at 850/im for dust 
radiating at Tj, k^^q is the dust emissivity at 850 fim (assumed 
here to be k^so - 0.02 cm^ g '), 'Wh2 is the molecular hydrogen 
mass, and ju is the mean molecular weight (in units of mu^). 

As shown in Appendix A, significant CO line flux can con- 
taminate the broadband continuum measurements; therefore, it 
is important to correct the continuum measurements for the 
presence of line emission. Having both line and continuum 
measurements for each of the sources in Orion allows for a di- 
rect measure of the level of this line contamination. At 850 fim, 
1 Jy within the JCMT 14 arcsecond beam is equivalent to a 
brightness temperature of 5 1 mK. However, the broad SCUBA 
850 fj.m bandpass, 30 GHz wide, dilutes the importance of indi- 
vidual spectral line features against the continuum. The equiv- 
alent line flux required to mimic a 1 Jy signal with SCUBA on 
the JCMT is 1540Kkms"'. For a typical peak gas tempera- 
ture of ~ 25 K, the Gaussian component of the CO line would 
need to be ~ 6kms ' wide to produce a ten percent, 0.1 Jy, 
error in the continuum measurement. Thus, for strong contin- 
uum sources, S > 2 Jy, line contamination should be minimal. 
Table|5]presents the CO line contamination for each of the six 
Orion sources. For completeness, the combined line contami- 
nation due to all other observed lines lying within the 850 fim 
bandpass is also shown. With the exception of the shock knot 
SKI, the total line contamination is always less than about ten 
percent. Note that SKI has a high peak brightness temperature 
for the CO 3-2 transition and has very extended Une wings, 
consistent with its identification as a hot shock front by Yu et 
al. (1997) on the basis of bright H2 emission. 
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After accounting for line contamination, it is straightfor- 
ward to compute the column density of H2, derived from the 
strength of the remaining continuum emission, and the frac- 
tional abundance of CO for each source. These values are also 
presented in Table |5] Within the PDR and the shock front, the 
abundance of gas-phase CO agrees with the expected value 
of 10""^ to within the accuracy of the computations (a factor 
of a few). However, toward each of the sub-millimeter proto- 
star regions the CO gas-phase abundance is smaller than ex- 
pected by about a factor of 5 - 10, further evidence that al- 
most all of the CO has been depleted onto grain surfaces in the 
dense envelopes. This result is in qualitative agreement with the 
more detailed modeling results presented by J0rgensen et al. 
(2002), where low-mass Class and pre-stellar sources were 
observed to have similarly depleted CO. Here, however, the 
outer envelope is warm, ~ 30 K and the bulk of the de- 
pletion likely occurs in an intermediate zone between the hot, 
centrally heated deep interior and the externally heated surface 
layer. The higher abundance of CO in MMS9 is likely due to 
the smaller envelope, and thus a fractionally smaller interme- 
diate zone. Detailed modeling of externally heated protostellar 
envelopes supports this hypothesis (J0rgensen et al. in prep.). 

4.3. Formaldehyde 

One of the advantages of H2CO is the ability to observe, in the 
same side-band of a single spectrum, two transitions which to- 
gether probe the temperature of the gas under a large variation 
in density conditions. This feature removes many of the obser- 
vational uncertainties inherent in sub-millimeter spectroscopy 
such as beam effects, weather conditions, and instrumental 
gain. The RADEX calculations for H2CO include energy lev- 
els up to y = 10 and use collisional rate coefficients by Green 
(1991). 

In the present study we observed three para-formaldehyde 
lines. Two of the lines were simultaneously observed in the 218 
GHz atmospheric window, 3o3 - 2()2 and 322 - 221. Because 
transitions between the different Kp ladders are radiatively for- 
bidden, the ratio of the line strengths from transitions in dif- 
ferent ladders depends only on the kinetic temperature of the 
gas, operating as a finely tuned thermometer over an extended 
density range (Mangum & Wootten 1993). A third line was ob- 
served at 362 GHz, 5o5 - 4o4, within the same Kp ladder as the 
3o3 - 2o2 transition. Within a given Kp ladder the excitation en- 
ergies do not vary dramatically and thus the ratio of these lines 
provides a useful density indicator at low to moderate densities 
(«H < 3 X 10^ cm""*) and a second temperature probe at higher 
densities. The first panel in Fig.|5lplots the expected line ratios 
against temperature and density. 

The physical parameters derived from the formaldehyde 
lines are presented in Table |6l and are plotted graphically for 
each source in Fig. [S] The spread of line ratios 7^33 and ki- 
netic temperature of the gas inferred from the 3o3 - 2o2 and 
322 - 221 transitions ranges from /?33 < 2, T > 300 K, at 
the location of PDRl, to R32 > 10, T < 30 K at the location 
of the sub-millimeter source MMS9. As expected, the inter- 
mediate mass star-forming regions NGC207 1 and FIR4 show 



moderately high temperatures T > 70 K, consistent with the 
observability of infrared dust continuum emission. Mangum & 
Wooten (1993) also observed FIR4 and NGC2071, obtaining 
temperatures of 85 K and 80 K and with a range of probable 
temperatures bounding the results of this study. The only in- 
frared source in OMC3, MMS7, shows no evidence for high 
temperatures. 

The strength of the 5o5 - 4o4 transition constrains the den- 
sity of the environment from which the emission is observed. 
A complication arises from the difference in the beam size at 
the telescope for the 3o3 - 2o2 and 5o5 - 4o4 lines. If the line 
emission fills the larger 218GHz beam, then there are no diffi- 
culties in taking the ratio of the brightness temperatures of the 
two lines. However, if the source of the formaldehyde emission 
is smaller than the 362GHz beam the ratio of the beam sizes 
should be taken into account when comparing the brightness 
temperatures to account for the effect of beam dilution. Source 
sizes intermediate between the 218GHz and 363GHz beam re- 
quire a varying normalization. Thus, the measured line ratio 
/?53, computed from the 5o5 - 4o4 and 3()3 - 2o2 lines, has an 
inherent uncertainty of ~ 3 due to the possible complications 
from intermediate beam dilution, with the value given in Table 
|6lstrictly an upper limit. 

For all sources the measured ratio /?53 is larger than 0.3, im- 
plying a high density, « > 10^ cm"^. The relatively evolved far 
infrared source FIR4 shows two kinematically distinct compo- 
nents, a narrow-line region (AV ~ 1.5 km s"') requiring an ef- 
fective density near 10^ cm"^ and broad-line component {AV ~ 
6kms"') requiring a much higher density. While MMS6 and 
NGC2071 also show distinct velocity components, the den- 
sity does not obviously vary between the components for these 
sources. Intriguingly, in all three cases the temperature deduced 
from both the narrow and broad-line components are quite sim- 
ilar 

With the temperature and density estimated, the column 
density of para-formaldehyde can be deduced from the inte- 
grated line strength of the 3o3 - 2o2 transition. The high temper- 
ature equilibrium ratio of the para to ortho forms (1:3) was used 
to convert this result into a total formaldehyde column density. 
The inferred abundance of formaldehyde was then measured 
against both the H2 column, deduced from the sub-millimeter 
dust continuum emission, and the measured CO column. The 
results are presented in Table|6l The two infrared sources show 
the largest abundance of formaldehyde measured against the 
CO. For these sources the abundance ratios lie in the range 
[H2CO]/[CO] = 2.5 - 10.0 X 10"^ The sub-millimeter sources 
MMS6 and MMS9 show a lower abundance ratio [H2CO]/[CO] 
~ 1 X 10"^'. Disentangling the importance of CO depletion in 
these data complicates any analysis unless the formaldehyde 
depletes in a similar fashion to CO. However, comparing the 
formaldehyde abundance to H2 produces a clear relation (see 
Table |6l and Fig. |8}. The hot PDR region and the infrared 
sources, NGC2071 and FIR4, have the largest formaldehyde 
abundance, approximately 5-10 times greater than observed for 
the sub-millimeter sources. At least half of the increase in the 
abundance of formaldehyde in the two infrared sources is due 
to the broad-line component. 
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In low-mass Class sources, typical [H2CO]/[CO] abun- 
dance ratios in the outer region are a few xlO"^, similar to those 
observed here, but jumps in the H2CO abundance by more than 
a factor of 100 have been found in the inner warm envelopes, 
consistent with ice evaporation (e.g., Ceccarelli et al. 2000, 
Schoier et al. 2002, Maret et al. in prep.). H2CO is expected 
to be a product of grain surface hydrogenation of CO and has 
been tentatively detected in interstellar ices with an abundance 
of a few percent with respect to H2O ice, corresponding to 
~ 10"^ with respect to H2 (Keane et al. 2001b). Thus, the trend 
observed here for the protostellar sources could be consistent 
with evaporation of H2CO-containing ices in the warm interior 
regions, with more energetic sources containing a larger frac- 
tional mass of envelope at high enough temperatures to evapo- 
rate such ices. For high-mass sources, H2CO shows a less clear 
relation with temperature (van der Tak et al. 2000b). 

4.4. Methanol 

Methanol emission lines have been detected in complex bands 
at 241 and 338 GHz, with excitation energies ranging upward 
from Tex ~30 K. The detection of many CH3OH lines in our 
spectra allows for the construction of a rotation diagram for 
each source in order to get a first indication of the tempera- 
ture and column density. Described in detail by Blake et al. 
(1987), this method assumes that the lines are optically thin 
and that the excitation can be characterized by a single tem- 
perature. The resulting rotation diagrams, however, are often 
misleading. Direct computation of the level populations using 
RADEX shows that the higher excitation lines require increas- 
ingly larger critical densities and thus the derived temperature 
from rotation diagram fitting is heavily dependent on the den- 
sity of the medium in which the methanol resides (see also the 
discussion by Bachiller et al. 1995). Figure |6l reveals the ef- 
fect on the rotation diagram produced by increasing the critical 
density of the medium while holding the kinetic temperature 
constant at 100 K and the total column density at A^chiOh = 
lO''* cm"^. Even though reasonable linear fits are produced, the 
derived temperature ranges from 25 K for = 10^ cm"^ to 
80 K for - 10^ cm"^. The derived column densities also 
overestimate the actual amount of methanol by factors of a few. 

For a protostellar source with both density and tempera- 
ture gradients a proper analysis becomes even more complex. 
Given that higher excitation lines tend to have higher critical 
densities, it is likely that these lines are dominated by emission 
solely from the deep interior of the envelope. Indeed, there is 
some evidence that the source size decreases as the required 
excitation of the line increases (Hatchell et al. 1998). The bulk 
properties of the envelope derived from methanol, however, are 
expected to be only slightly weighted toward the highest den- 
sity and temperature regions due to the small angular extent of 
these extreme locations. Complicating even a detailed analy- 
sis, there is evidence for an enhanced abundance of methanol 
at high temperatures, T > 90 K, in protostellar envelopes (van 
der Tak et al. 2002b, Schoier et al. 2002). 

Despite the above concerns, one might hope that methanol 
may be used both as a temperature and density indicator by 



treating all of the lines through a proper statistical equilibrium 
calculation. The RADEX calculations for methanol include en- 
ergy levels up to 260 K and use collisional rate coefficients by 
Lees & Haque (1974) and Walmsley (private communication). 
For this paper, all the measured methanol lines were compared 
with RADEX computations to produce a chi-squared fit as a 
function of temperature and density. The results are presented 
graphically in Fig. and are tabulated in Table For FIR4, 
the methanol measurements can be separated into narrow and 
broad line components. The narrow line region requires very 
high density and temperature n > 10^ cm"^, T ~ 7QK, while 
the broad line component requires n ~ 10'^ cm"^, T ~ 80 K. 
The broad line region appears to be beam diluted, at least with 
respect to the larger 241 GHz beam, requiring a normalization 
factor between the two sets of measurements to produce rea- 
sonable chi-squared fits. These results are not consistent with 
the formaldehyde measurements in that higher densities are re- 
quired to account for the methanol observations. Given a range 
of densities within the surroundings of FIR4, however, the 
methanol observations would be most sensitive to the highest 
density gas. In contrast, the other infrared source, NGC2071, is 
best fit by a moderate density, n ~ 5 x 10^ cm"-', and a mod- 
erate temperature, T ~ 60 K, consistent with the formaldehyde 
measurements. 

The two sub-millimeter sources show some variation in the 
infeiTed temperature. MMS6 requires a moderate temperature 
r ~ 50 K while MMS9 requires a low temperature T < 30 K. 
For both sources, the methanol appears to probe a dense in- 
ner region, ~ 10^ cm"^ compared with the formaldehyde de- 
rived envelope densities ~ 10^ cm"^. This is likely due to the 
higher critical density required to excite the methanol transi- 
tions. Across all the sources there is a reasonable correlation 
between the methanol and formaldehyde temperatures with the 
sub-millimeter sources being the coldest, the infrared sources 
being moderately warm, and the PDR having the highest tem- 
peratures. 

The derived abundance of methanol is similar to that of 
formaldehyde. The warm infrared sources have the highest 
abundance with respect to CO. However, when compared with 
the total column of hydrogen deduced from the sub-millimeter 
continuum measurements (Table|5}, the PDR region also shows 
an enhanced abundance of methanol. Considering that the 
methanol observations are predominantly probing higher den- 
sity regions within each envelope, it is likely that these abun- 
dance estimates are strictly lower limits. Differences in the 
abundance derivations across sources may be affected by the 
extent of the warm interior within each envelope. 

Like H2CO, CH3OH is expected to be a product of the hy- 
drogenation of CO on grain surfaces (Tielens & Charnley 1997) 
and is detected in grain mantles, albeit with significant varia- 
tions between different sources (Dartois et al. 1999). CH3OH 
abundance jumps of at least a factor of 100 have been found in 
various high-mass objects (van der Tak et al. 2000b), whereas 
the well-studied low-mass Class source IRAS 16293 -2422 
shows a jump by nearly three orders of magnitude in the inner 
warm envelope (Schoier et al. 2002). Such high CH3OH abun- 
dances can only be produced by evaporation from grain man- 
tles. Enhanced CH3OH abundances have also been detected in 
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outflows at positions offset from the protostar, e.g., in L1157 
(Bachifler & Perez-GutieiTez 1997), where the CH3OH is re- 
leased from the grain mantles by shocks rather than thermal 
evaporation. The general trend of an increased CH3OH abun- 
dance with temperature observed for our sources is consis- 
tent with this picture. The presence of both narrow and broad 
CH3OH line profiles suggests that both evaporation mecha- 
nisms play a role in some objects (e.g., NGC 2071). 

4.5. Nitrogen-bearing Species 
4.5.1. HCN and HNC 

Identical lines of both HCN and HNC are observed simulta- 
neously with a single spectrograph setting at the JCMT and 
thus, assuming the collisional cross sections for the various ex- 
citations are identical across species, the relative abundance of 
the two molecules can be simply determined if the transitions 
are optically thin. To test the optical depth conditions, the ra- 
tios of line strengths between HCN and H'^CN, and HNC and 
HN'-'C were computed. In all cases the HCN line was found 
to be moderately to strongly optically thick (thcn ~ 3 - 6) 
however, only the most energetic sources FIR4 and NGC207 1 
were found to have optically thick HNC lines (thnc ~ 2 - 3). 
The [HCN]/[HNC] ratio tabulated in Table E] was determined 
by accounting for these optical depth efi'ects. 

For the four protostellar regions observed, the integrated 
[HCN]/[HNC] abundance ratio ranges from 5-10. Using the 
chemical models by Schilke et al. (1992), this would imply 
low temperatures T < 30K and high densities n > 10^ cm"^, 
or high temperatures T > 50 K and moderate to low densities 
n < 10*' cm"^. For PDRl, the ratio rises to 70, requiring both 
high densities n > 10^ cm"^ and high temperatures T > 50K 
(Schilke et al. 1992). It is worth noting that the main obser- 
vational difference between the least energetic sub-millimeter 
sources, MMS6 and MMS9, and the most energetic sources, 
FIR4 and NGC2071, is in the peak brightness of the lines. 
Thus, it would appear that the conditions probed by the HCN 
and HNC molecules in all four sources are nearly identical 
but that the solid angle covered in the most energetic sources 
is somewhat larger. RADEX calculations support this hypoth- 
esis. For the determined physical conditions of the optically 
thick HCN lines the excitation temperature is expected to be 
relatively large Tex > 10 K. The least energetic sources have 
observed peak intensities of only Tpeak ~ 1 K, while the in- 
frared sources show Tpeak ~ 15 K. Care should be taken in us- 
ing chemistry to constrain temperature, however, since the low- 
mass sources observed by J0rgensen et al. (2003; in prep.) have 
high HCN/HNC abundance ratios at low temperatures, which 
are not readily explained by the Schilke et al. (1992) models. 

RADEX calculations for HCN, including energy levels 
up to y - 12 and using collisional coefficients by Green & 
Thaddeus (1974), have been performed to derive HCN abun- 
dances using the estimated physical conditions derived from 
formaldehyde. The resulting abundances with respect to CO are 
~5xlO"^ for the sub-millimeter sources and PDRl; however, 
the production of HCN appears to be significantly enhanced in 
the far infrared sources FIR4 and NGC2071. The HCN abun- 



dances with respect to H2 vary from 0.6 x 10"'" to 40 x 10"'" 
(Table |8|i. These latter abundances show a clear trend with in- 
creasing source energetics. HCN has not been detected in inter- 
stellar ices due to a lack of clear spectroscopic signatures, but is 
expected to freeze-out in cold regions like most molecules. The 
HCN abundance is known to be enhanced by orders of magni- 
tude in the warm inner regions of massive protostars (Lahuis & 
van Dishoeck 2000; Boonman et al. 2001), but the high temper- 
ature routes to HCN do not become effective until T > 200 K 
(e.g.. Doty et al. 2002). It remains to be determined whether the 
intermediate mass sources studied here have sufficient amounts 
of gas at such high temperatures to enhance HCN significantly. 
The trends observed here are more likely due to general evap- 
oration of ices with temperatures between 20 and 100 K, en- 
hancing the overall gas-phase abundances. 

Detailed consideration of the HCN line shows that it is usu- 
ally wider than the HNC line. This requires that in the broad 
line region the ratio of the two molecules is much higher than 
the ratio of the total integrated intensities, implying higher tem- 
peratures and/or higher densities if chemical equilibrium is as- 
sumed (Schilke et al. 1992). Another possible explanation for 
the general lack of broad line HNC may be found in the shock 
model solution of Schilke et al. (1992) where HNC is pref- 
erentially destroyed within the interaction region, reforming 
downstream from the shock front. In particular, the broad HCN 
emission, coupled with little or no HNC emission, found in 
the shock knot SKI is likely best explained in this manner. 
Interferometer data of the low-mass source IRAS 16293 -2422 
also show that HNC avoids the warm inner region, in contrast 
with HCN (Schoier et al. 2003; in prep.). 

4.5.2. CN 

Observations of the CN hyperfine structure at 340 GHz pro- 
vides a set of constraints for the modeling of CN (see Table|3}- 
In the optically thin limit, the ratio between the strongest and 
weakest hyperfine components should be ~ 15; however, this 
is only observed in the spectrum of PDRl . In all other regions, 
except possibly MMS6, the main CN component is optically 
thick, (tcn ~ 3 - 6), while the weaker components remain rea- 
sonably optically thin. Assuming that the line flux is predom- 
inantly from the outer envelope surrounding the sources, the 
column of CN derived from the optically thin hyperfine lines is 
presented in Table|8] The RADEX calculations for CN include 
energy levels up to 7 - 15. As found for the abundance of 
HCN, only the coldest, least energetic source, MMS6 has a dra- 
matically low abundance of CN. However, a clear correlation 
is found between the ratio [HCN]/[HNC] and [HCN]/[CN]. 
Assuming that all CN is locked up in these three molecules, 
the total abundance of CN can be readily determined. Table |8] 
displays the results of this calculation. With the exception of 
the sub-millimeter source (MMS6) the derived total abundance 
of CN is much more constant than any of the individual molec- 
ular species. 

Results of systematic studies of CN in low- or high-mass 
sources are not yet available, but like HNC, CN seems to avoid 
the inner warm region in the low-mass protostar IRAS 16293 
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-2422. The radical is known to be abundant in the outer regions 
of PDRs (e.g., Fuente et al. 1998), where it is produced by rapid 
gas-phase chemistry of CH and C2 with N. Alternatively, it can 
be a photodissociation product of HCN. Photodissociation of 
CN itself requires photons with energies higher than 12 eV. 
The observed lack of variation of the CN abundance in the pro- 
tostellar sources observed here implies that they are exposed to 
similar radiation fields. As mentioned earlier, the Orion region 
is bathed in a significantly increased external radiation field due 
to the close proximity of massive stars. The much lower abun- 
dance within the massive envelope of MMS6, however, is not 
easily explained via this hypothesis unless the structure of the 
envelope is such that the interstellar radiation field is unable 
to penetrate through the surface layer and into the significant 
interior. 

4.6. Sulphur-bearing Species 

4.6.1. CS 

Both CS and the rarer C^"*S isotope were observed in the 7 = 5- 
4 transition. In all cases the line ratio was less than the iso- 
topic ratio of 22 implying moderately optically thick conditions 
within the CS line but optically thin conditions for the isotope. 
Taking the formaldehyde physical conditions for the tempera- 
ture and the density, the column density of CS was determined 
using RADEX (Table|9}. The RADEX calculations for CS in- 
clude energy levels up to 7 = 12. The abundance with respect to 
H2 trends with the source energetics, rising from a few 10 " to 
a few 10 in the protostellar sources. This is a stronger trend 
than seen for HCN suggesting that the more energetic sources 
are providing more heating to their outer envelope. The abun- 
dance rises sharply to a few 10"^, in the warm PDRl source. 

For four of the sources the J = 7-6 transition of C^'^^S was 
also detected. In all observed cases the ratio of the J = 7-6 
versus J = 5^ lines implies a warm, (T > 50 K) and dense 
(n > 10^) environment (Table |9}. The C^'*S-derived temper- 
ature is higher than that derived from H2CO and more likely 
refers to the inner warm and dense region rather than the outer 
envelope. If more transitions are observed, a detailed density 
structure for each of the sources could be determined. 

4.6.2. SO 

The only observed SO transition was unfortunately blended 
with a methanol line and thus difficult to accurately measure. 
Despite this complication, three sources, spanning a range of 
energetics and conditions, have unambiguous detections for 
SO (MMS6, NGC2071, and PDRl). Assuming the formalde- 
hyde derived properties for the physical conditions of the gas, 
RADEX was used to determine the required column density of 
SO for each of the sources. The RADEX calculations for SO 
include energy levels up to 580 K. The results are presented in 
Table|5] along with the abundance of SO in each source and the 
SO to CS abundance ratio. 

Observations of more SO lines and other sulphur-bearing 
species, such as SO2 and H2S, are needed to study possible 
differences in the sulphur-chemistry in these sources. 



5. Discussion and Conclusion 

This molecular line study was undertaken to determine if mor- 
phological clues and qualitative indicators were observable 
across a range of environmental conditions, from pre-stellar 
and young protostellar envelopes Lboi ~ 1 - IOOL0 through 
infrared bright energetic sources Lboi ~ 4OOL0, a PDR knot 
and a shock front. Several trends are apparent, especially in the 
derived abundances of many molecular species. As well, there 
are a number of spectral hints that protostellar sources reside 
within the sub-millimeter clumps, excluding the PDR and the 
shock knot. 

5.1. Abundances 

Despite the inherent danger in assuming a single temperature 
and density for the environments of each source, general trends 
in abundance are observed. It is worth noting that detailed mod- 
eling (van der Tak et al. 2000a, Schoier et al. 2002, Doty et al. 
2002, J0rgensen et al. 2002) provides a much more accurate 
determination of relative abundances, especially for molecules 
which are excited only in parts of the envelope. Also, the pres- 
ence of abundance gradients or 'jumps' can be established for 
some molecules. Such models require a determination of the 
density and temperature profiles of the sources from the dust 
continuum emission, which has not yet been completed for 
this study. Schoier et al. (2002) show that the inferred abun- 
dances using a constant temperature and density do not differ 
by more than a factor of a few from the detailed analysis as 
long as the adopted conditions are appropriate for the partic- 
ular molecule or line, which should be the case for this work. 
However, detailed modeling requires significant input as to the 
density and temperature profile of the source, as well as relying 
on additional assumptions such as the dust emissivity profile. 
Additionally, the Orion sources are bathed in a strong external 
radiation field (Li et al. 2003) requiring careful consideration of 
the exterior conditions of each source envelope where the dust 
and gas temperatures reach T ~ 30K. This study is concerned 
primarily with the constraints on physical conditions provided 
by the variations in the molecular tracers without resorting to 
detailed modeling of individual sources. 

Some strong trends are observed across the source list. The 
peak brightness and integrated line strength of both CO and 
'^CO 3-2 lines follow the energetics and warm dust tempera- 
tures of the envelopes. Despite the observed depletion of CO in 
the protostellar source envelopes (typically the abundance is a 
factor of 10 lower than the mean molecular cloud value), in all 
cases the objects are visible in each of the observed CO lines. 

The formaldehyde derived temperatures, T(H2CO), do not 
match the envelope dust temperatures, measured in the sub- 
millimeter However, for the sub-millimeter sources MMS6, 
MMS7, and MMS9 T(H2CO) is consistent with the gas tem- 
perature in the outer envelope as seen in the '^CO line. For the 
infrai-ed sources, FIR4 and NGC2071, T(H2CO) requires in- 
ternal heating through a significant fraction of each envelope. 
Detailed modeling of other sources provides evidence that the 
formaldehyde abundance increases in the warm, dense interi- 
ors of protostellar envelopes, and that grain mantle evaporation 
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may be important in producing the enhancement. However, the 
relatively low lying lines observed in this study are more accu- 
rate tracers of the extended envelope conditions, where most of 
the mass resides (J0rgensen et al. in prep.). While more uncer- 
tain, the methanol derived temperatures follow the same trend 
as the formaldehyde temperatures. For the remainder of the dis- 
cussion we adopt the formaldehyde conditions as representa- 
tive of the bulk conditions within each source. 

Despite uncertainties of order unity, the abundance of both 
formaldehyde and methanol correlate well against the source 
energetics (Fig.|^. The changing abundance may reflect jump 
conditions in formaldehyde and methanol abundance with dust 
temperature in the inner warm part. Ignoring the broad (out- 
flow) methanol component in NGC207 1 , the abundance of both 
formaldehyde and methanol may be plotted as a step func- 
tion with cold sources, T < 50 K, having a low abundance 
and warm sources, T > 50 K, having about five times higher 
abundance. Although less observational data from the study 
are available, the CS abundance also appears to correlate with 
source energetics. 

The abundance of HCN is also sensitive to source con- 
ditions; however, this is mostly determined by extreme cases 
MMS6 and PDRl. In contrast to the other observed species, 
little variation of abundance was found for HCN among the 
energetic sources, F1R4 and NGC2071, and the weak sub- 
millimeter source MMS9. Considering the three dominant CN- 
bearing species (HCN, HNC, and CN), the total abundance of 
these three molecules is quite constant across a wide range of 
conditions, varying by less than thirty percent between MMS9, 
F1R4, and NGC2071 and only varying by a factor of two 
when the PDR region is included. Only the cold, dense sub- 
millimeter source MMS6 appears to have a severely depleted 
abundance of CN-bearing species. The enhanced external radi- 
ation field in Orion may be responsible for this apparent equi- 
librium. 

The general observed trend of increased abundance with 
source energetics is consistent with a scenario in which freeze- 
out of molecules occurs in the cooler intermediate zone within 
the envelope and evaporation of ices occurs in the warmer inte- 
rior and possibly exterior regions. Similar trends have been ob- 
served for samples of low-mass (J0rgensen et al. 2002, 2003) 
and high-mass (van der Tak et al. 2000a,b; Boonman et al. 
2003) sources without external heating, although not for all 
molecules. The only other intermediate-mass source that has 
been studied in some detail chemically is AFGL 490 (Schreyer 
et al. 2002) which has a luminosity of lO-' L©. This source has 
a large envelope, with abundances comparable to those of the 
warmer sources observed here. Strong solid-state features of 
various species indicative of freeze-out are also detected. 

None of the observed sources shows the characteristic 
crowded line spectra of a 'hot core', such as found for high- 
mass protostars like G34.3+0.15 (Macdonald et al. 1996), 
G327.3-0.6 (Gibb et al. 2000) or W3(H20) (Helmich & van 
Dishoeck 1997). Deeper integrations are needed to determine 
whether complex organic molecules like CH3OCH3 are present 
in the warmer sources in our sample, e.g. FIR4 or NGC 2071. 
These molecules are expected to be produced by gas-phase 
chemistry between evaporated ices, but have so far only been 



observed in high-mass sources. The different dynamical time 
scales of the hot core gas in low- versus high-mass objects com- 
pared with the chemical timescales of ~ 10"* yr may prevent 
formation of such second generation species. 

The PDRl position differs chemically from the other 
sources by having the highest gas-phase abundances overall 
and the largest [HCN]/[HNC] ratio. PDRs also generally have 
high abundances of radicals such as CN and C2H and ions like 
CO^ (Jansen et al. 1995). Indeed, a study of a set of more 
evolved intermediate-mass sources by Fuente et al. (1998) has 
identified CN as a particularly good diagnostic of the ultravio- 
let radiation. In our sources, this trend is not so evident, likely 
because all sources are located in the Orion region which is 
bathed in intense radiation from nearby O and B stars. 

Few molecules are observed in the SK1-OMC3 shock and 
only CO shows truly broad line wings at that position. Sulphur- 
containing species like H2S, SO and SO2 are predicted to be 
abundant in shock models (e.g., Leen & Graff 1988) and are 
seen to be enhanced in the Orion-KL plateau gas with broad 
line wings, but are not prominent here; however, the source is 
quite weak and the lines may be below our detection limits. In 
NGC 2071, molecules like CH3OH are present in the outflow 
but this results from grain evaporation in the shock rather than 
high-temperature chemistry. Chernin, Mason, & Fuller (1994) 
have observed broad SO emission in NGC 2071 with some 
abundance enhancement. Deeper searches for shock diagnos- 
tics and accurate determinations of abundances at shocked and 
non-shocked positions in these and other sources are needed. 

5.2. Protostellar Source Diagnostics 

A large number of sub-millimeter continuum surveys have now 
been completed in star-forming molecular clouds, and hun- 
dreds of new sub-millimeter bright envelopes have been enu- 
merated. However, an outstanding question is which fraction 
of these objects surround protostellar sources and which are 
pre-stellar Within the present study there were no clear pre- 
stellar objects although MMS6 appear somewhat ambiguously 
defined as a Class source. Despite this lack of an obvious pre- 
stellar candidate, it is worth considering if any of the observed 
molecular signatures require an embedded heating source. 

For many of the sub-millimeter objects found within sur- 
veys, '^CO observations have provided no indication of a co- 
incident CO peak (Mitchell et al. 2001). Thus, the clear mea- 
surement of CO isotopes in this sample, despite the depletion, 
may provide a clue to which objects contain embedded sources, 
perhaps by warming and evaporating CO back into the envi- 
ronment from which it had frozen out during an earlier epoch. 
Alternatively, the presence of a strong external radiation field 
may bias the results in Orion. The broad CO line wings provide 
clear evidence of enhanced kinematics within these condensa- 
tions and, along with the ubiquitous outflows in star-forming 
regions, act as a sign post for embedded sources. 

The appearance of molecular lines with relatively high 
excitation temperatures, and preferably also with high criti- 
cal densities such that the warm region might be inferred to 
be deeper within the envelope, should provide a signature to 
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Class and later sources. Thus, the methanol lines are likely 
the strongest indicator of a warm, dense region within the 
envelope. For all sources observed in this study, except the 
known outflow source MMS9, the required internal tempera- 
ture deduced from the methanol lines is > 50 K, a temperature 
unattainable without an energetic internal source. The CS 7- 
6 line also provides evidence of a warm, dense interior region 
in all sources except MMS9. It is interesting to note that all 
the protostellar sources considered in this paper require high 
densities in their interior Standard models for star formation 
expect a power-law density distribution toward the clump cen- 
ter during the collapse phase (Shu 1977; Henriksen, Andre, & 
Bontemps 1997) The high densities measured here may be a 
result of the stage of evolution of the individual clumps. 

Formaldehyde observations also provide important evi- 
dence for a large warm envelope. While the low-lying lines 
observed in this study are predominantly excited in the mas- 
sive outer envelope, the derived temperature places constraints 
on the required heating source. In particular, envelope temper- 
atures above 20 K are difficult to reconcile without either a 
warming source inside or an external heating source. All the 
sources observed in this study have two components, narrow 
and broad, associated with the formaldehyde measurements. 
However, this is not always the case for sub-millimeter bright 
regions (Mitchell et al. 2001; Tothill, private communication). 
It is possible that the broad component of the formaldehyde, 
often implying a much denser zone, is tracing an inner region 
within the envelope which is undergoing collapse or has be- 
come much more turbulent. 

Not all heating is due to protostellar sources and it is there- 
fore important to distinguish external heating (such as in a 
PDR) from internal heating. Protostellar sources should have 
central condensations which show up both in sub-millimeter 
dust continuum maps and maps of molecular lines with high 
critical densities (e.g., CS 7-6). Chemically, PDRs are best rec- 
ognized by a high [HCN]/[HNC] ratio and by high abundances 
of radicals such as CN and C2H and ions like CO^ (see above 
discussion). 

5.3. Line Contamination 

One of the original motivations for this molecular line study 
was to determine the importance of Une contamination within 
the broad SCUBA 850yum passband. While theoretical calcu- 
lations (Appendix A) show that the influence of CO contami- 
nation can become exceedingly large in particular situations, it 
is also clear that typical conditions within molecular clouds are 
not so extreme. The observational evidence presented (Table 
|5} does show tiiat the contamination, while typically less than 
10%, occasionally dominates the continuum flux. However, 
this only occurs in regions with warmer molecular gas temper- 
atures and large velocity gradients which allow for enhanced 
integrated line strengths in CO and its isotopes. Such con- 
ditions occur most often at isolated shock fronts within the 
cloud, for example in the knots associated with protostellar jets. 
Contamination may also arise from lines of other molecules, 
especially around energetic sources. The best observed exam- 



ple for this is Orion KL, for which a forest of lines, primarily 
from SO and SO2, produce between 28% and 50% line con- 
tamination at 850/im (Serabyn & Weisstein 1995; Groesbeck 
1994). Consideration of Table|5]shows that in our sample the 
line contamination is never dominated by lines other than CO; 
however, for the most energetic sources line emission from 
HCN, HNC, CN, and methanol provides a significant fraction 
(> 40 %) of the total line contamination. 
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Name R.A}"^ DccJ"' Tj M,,,,. L Object 



(J2000) (J2000) {K) (Mo) (Lq) 



MMS6 - OMC3 


5 : 


: 35 


: 23.5 


-05 : 


: 01 


: 32 


15 


36 


1.2-60 


ClassO 


MMS9 - 0MC3 


5 : 


: 35 


:26.2 


-05 : 


: 05 


:46 


20 


10 


0.6 - 90 


ClassO 


MMS7 - 0MC3 


5 : 


: 35 : 


: 26.4 


-05 : 


: 03 : 


; 59 


26-50 


8 


76 


ClassO 


FIR4 - 0MC2 


5 : 


: 35 : 


: 26.7 


-05 : 


: 09 : 


; 59 


15-40 


34 


400 


ClassO/I 


NGC2071IR 


5 : 


:47; 


: 04.4 


+00: 


: 21 : 


: 49 


20-50 


30 


500 


ClassO/I 


PDRl - OMCl 


5 : 


: 35 


; 25.3 


-05 : 


: 24 


: 34 


60 






PDR 


SKI - 0MC3 


5 : 


: 35 ; 


: 17.0 


-05 : 


:06; 


: 03 


10 






Shock 



Positions are taken from the sub-millimter maps of the ISF 
(Johnstone & Bally 1999) and Orion B North (Johnstone et al. 
2001). 
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Table 2. Calendar of Observations and Conditions 



Date 


Receiver 


^upper 


^lower 








Lines 






(GHz) 


(GHz) 


(km/s) 


(K) 




1999 January 


RxA3 


218.35 


226.48 


0.2 


400- 


700 


HzCO.CHjOH 


1999 September 


RxB3 


330.53 


338.53 


0.5 


1200- 


1700 


"CO.CHsOH 


1999 September 


RxB3 


338.53 


346.53 


0.5 


1200- 


1700 


CH30H,SO 


1999 September 


RxB3 


354.40 


362.40 


0.5 


1200- 


1700 


HCN,HNC 


2000 December 


RxA3 


233.00 


241.00 


0.4 


350- 


400 


C'^S 


2000 December 


RxA3 


233.75 


241.75 


0.4 


350- 


400 


CH30H 


2000 December 


RxA3 


236.94 


244.94 


0.4 


375- 


385 


cs 


2000 December 


RxB3 


337.25 


345.25 


0.5 


530 - 550, 


710-740 




2000 December 


RxB3 


340.45 


348.45 


0.5 


700- 


800 


CN, HN"C 


2000 December 


RxB3 


345.67 


353.67 


0.5 


490- 


■510 


CO 


2000 December 


RxB3 


354.60 


362.60 


0.5 


590- 


620 


HCN.HNCHjCO 
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Table 3. Observed Spectra and Gaussian Fit Parameters 



Source Species 


V 


J 


W 


AV 


Z'jVIB 




(GHz) 




(Kkms-') 


(kms-i) 


(K) 


MMS6-OMC3 CO 


345.796 


3-2 


93.00 






"CO 


330.588 


3-2 


4.22 


0.57 


7.00 








40.70 


2.24 


17.1 


c"o 


337.061 


3-2 


4.78 


1.08 


4.15 


H2C0 


218.222 


303-202 


2.81 


0.88 


3.00 








2.52 


2.40 


0.99 


H2C0 


218.475 


3 -2 

22 21 


0.53 


1.73 


0.29 


H2C0 


362.736 


-'OS '04 


1.80 


0.91 


1.85 








1.76 


2.92 


0.57 


CH30H 


218.440 


42E-3iE 


0.44 


1.41 


0.29 


CH30H 


241.700 


5nE-4nE 


0.18 


0.82 


0.21 


CH30H 


241.767 


5_|E-4_iE 


0.95 


1.63 


0.55 


CH30H 


241.791 


5„A+-4nA+ 


1.27 


1.65 


0.72 


CH30H 


241.904 


5±2E-4^2E 


0.34 


2.80 


0.11 


CH30H 


338.345 


7_iE-6_iE 


0.84 


1.56 


0.50 


CH30H 


338.404 


7fiE-6fiE 


0.34 


0.99 


0.32 




338.409 


7nA+-6nA+ 


0.96 


3.45 


0.26 


CH30H 


338.583 


73E-63E 


0.51 


4.81 


0.10 


CH30H 


338.615 


7iE-6iE 


0.46 


2.05 


0.21 


CH30H 


338.722 


7±2E-6±2E 


0.45 


1.69 


0.25 


CN 


340.248 


3-2 


2.95 


1.41 


1.97 


HNC 


362.630 


4-3 


4.56 


1.06 


4.04 


HCN 


354.505 


4-3 


9.65 


3.81 


2.38 


H"CN 


345.340 


4-3 


0.49 


2.05 


0.23 


cs 


244.936 


5-4 


2.31 


0.93 


2.34 








2.23 


2.85 


0.74 


C^^'S 


241.016 


5-4 


0.65 


1.58 


0.39 


c^-'s 


337.396 


7-6 


0.54 


1.19 


0.43 


so 


346.528 


8,9-7,8 


1.28 


1.23 


0.98 
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Table 3. Observed Spectra and Gaussian fit parameters (cont'd) 



Source 


Species 


V 


J 


W 


AV 








(GHz) 




(Kkms-') 


(kms-') 


(K) 


MMS7-OMC3 


H2CO 


218.222 


3o3-2o2 


2.03 


1.38 


1.39 




H2CO 


218.475 




0.10 


0.55 


0.18 


MMS9-OMC3 


CO 


345.796 


3-2 


225.00 








"CO 


330.588 


3-2 


59.10 


2.40 


23.2 




c"o 


337.061 


3-2 


2.88 


1.29 


2.10 




H2C0 


218.222 




1.42 


0.67 


1.99 










1.65 


2.11 


0.73 




H2C0 


218.475 


J 22 -^21 


0.28 


0.74 


0.36 




H2C0 


362.736 


5o5-4o4 


0.94 


1.12 


0.78 




CH30H 


218.440 


42E-3iE 


0.23 


0.98 


0.22 




CH30H 


241.700 


5iiE-4nE 


0.14 


0.82 


0.16 




CH30H 


241.767 


5_iE-4_iE 


0.51 


1.21 


0.40 




CH30H 


241.791 


5„A+-4nA+ 


0.66 


0.86 


0.72 




CH30H 


338.345 


7_iE-6.iE 


0.30 


0.55 


0.50 




CN 


340.248 


3-2 


4.65 


1.78 


2.46 




HNC 


362.630 


4-3 


2.97 


0.98 


2.85 




HCN 


354.505 


4-3 


1.05 


0.82 


1.21 










3.57 


3.40 


0.98 




H"CN 


345.340 


4-3 


0.41 


5.23 


0.07 




CS 


244.936 


5-4 


2.73 


1.50 


1.71 
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Source Species 


V 


J 


W 




7'mb 








(^Jv Km s } 


\KSW S ) 




FIR4-OMC2 CO 


345.796 


3-2 


442.00 






"CO 


330.588 


3-2 


76.20 


2.57 


27.8 








24.70 


12.5 


1.85 


C"0 


337.061 


3-2 


1.27 


0.90 


1.32 








2.70 


1.65 


1.54 


H2CO 


218.222 


3o3-2o2 


9.82 


1.49 


6.21 








15.70 


6.28 


2.35 


H2CO 


218.475 


322-221 


2.96 


1.50 


1.85 








4.49 


5.81 


0.73 


H2CO 


362.736 


5 -A 


5.24 


1.20 


4.09 








29.3 


6.02 


4.50 


CH3OH 


218.440 


42E-3iE 


3.43 


0.98 


3.29 








7.49 


5.10 


1.38 


CH3OH 


241.700 


5nE-4nE 


2.39 


1.36 


1.65 








5.36 


5.33 


0.95 


CH3OH 


241.767 


5_iE-4-.iE 


4.73 


1.29 


3.44 








9.33 


5.50 


1.60 


CH3OH 


241.791 


5„A+-4nA+ 


5.93 


1.36 


4. 11 








10.10 


5.25 


1.81 


CH3OH 


241.904 


5+2E-4+,E 


2.35 


1.54 


1.43 








6.65 


5.68 


1.10 


CH3OH 


338.345 


7_iE-6_iE 


5.86 


1.34 


4.12 








21.70 


5.28 


3.85 


CH3OH 


338.404 


76E-66E 


5.31 


1.21 


4.14 








25.20 


5.30 


4.46 


CH3OH 


338.513 


74A±-64A=' 


1.58 


0.57 


2.59 








13.30 


7.37 


1.69 


CH3OH 


338.542 


73A*-64A* 


5.37 


2.29 


2.20 








11.30 


3.54 


3.00 


CH3OH 


338.583 


73E-63E 


1.44 


3.03 


0.45 








1.96 


5.70 


0.32 


CH3OH 


338.615 


7iE-6iE 


2.24 


1.26 


1.67 








10.20 


5.68 


1.69 


CH3OH 


338.640 


72A+-62OA+ 


9.44 


4.19 


2.12 


CH3OH 


338.722 


74.9E-6j-9E 


11.60 


2.70 


4.03 








29.20 


7.23 


3.79 


CN 


340.248 


3-2 


7.48 


1.88 


3.73 








15.50 


8.68 


1.67 


HNC 


362.630 


4-3 


9.63 


1.62 


5.60 








11.50 


6.83 


1.58 


HCN 


354.505 


4-3 


163.00 


10.2 


15.0 


h"cn 


345 340 


4-3 


8.66 


8.43 


0.97 


C^^S 


241.016 


5-4 


0.72 


1.40 


0.48 








1.31 


5.18 


0.24 


C^^S 


337.396 


7-6 


2.64 


7.39 


0.34 


so 


346.528 


8,9-7,8 


<9.00 


6.00 


0.97 








6.70 


6.00 


1.10 
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Table 3. Observed Spectra and Gaussian fit parameters (cont'd) 



Source Species 


V 


J 


W 


AV 


7'mb 




(GHz) 




(Kkms-') 


(kms"') 


(K) 


NGC2071IR CO 


345.796 


3-2 


770.00 






C"0 


337.061 


3-2 


2.96 


1.73 


1.60 








4.25 


7.31 


0.55 


H2CO 


218.222 


3o3-2o2 


12.00 


2.21 


5.10 








14.00 


9.12 


1.45 


H2CO 


218.475 




3.75 


2.95 


1.19 








3.67 


16.3 


0.21 


H2CO 


362.736 


-^05 '04 


10.90 


2.07 


4.93 








10.70 


7.50 


1.34 


CH3OH 


218.440 


42E-3iE 


3.12 


2.30 


1.28 


CH3OH 


241.700 


5nE-4nE 


2.39 


1.36 


1.65 








5.36 


5.33 


0.94 


CH3OH 


241.767 


5-iE-4_iE 


6.51 


2.38 


2.56 








6.18 


12.1 


0.48 


CH3OH 


241.791 


5oA^-4oA^ 


4.93 


1.80 


2.58 








8.12 


7.25 


1.05 


CH3OH 


241.879 


5iE-4iE 


1.45 


1.88 


0.73 


CH3OH 


241.888 


52A+-42A+ 


0.74 


2.24 


0.31 


CH3OH 


241.904 


54.9E-4^.9E 


1.09 


1.52 


0.67 








2.98 


4.66 


0.60 


CN 


340.248 


3-2 


14.50 


2.61 


5.20 








14.50 


8.63 


1.58 


HNC 


362.630 


4-3 


12.30 


1.92 


6.00 








10.70 


6.44 


1.56 


HN"C 


345.340 


4-3 


1.30 


3.89 


0.31 


HCN 


354.505 


4-3 


52.10 


3.95 


12.4 








31.20 


14.3 


2.05 


H"CN 


345.340 


4-3 


5.51 


10.1 


0.51 


CS 


244.936 


5-4 


10.30 


1.87 


5.17 








25.50 


6.29 


3.81 


C^S 


241.016 


5-4 


0.59 


0.85 


0.66 








3.86 


3.45 


1.05 


C'^^S 


337.396 


7-6 


0.53 


3.70 


0.13 








3.57 


3.74 


0.89 
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Source 


Species 


V 

(GHz) 


J 


W 

(Kkms-') 


AV 
(kms"') 


(K) 


PDRl-OMCl 


CO 


345.796 


3-2 


535.00 








"CO 


330.588 


3-2 


129.00 


1.39 


87.0 










90.10 


1.28 


66.0 




c"o 


337.061 


3-2 


10.20 


1.38 


6.94 




H2C0 


218.222 


3o3 2o2 


0.38 


0.49 


0.72 










5.11 


1.88 


2.55 




H2C0 


218.475 


322-221 


0.23 


0.49 


0.43 










1.76 


1.32 


1.25 




H2C0 


362.736 


5 -A 

•^05 '04 


3.80 


1.56 


2.29 




CH30H 


241.700 


5oE-4oE 


0.33 


1.83 


0.17 




CH30H 


241.767 


5_|E-4_,E 


0.89 


1.87 


0.45 




CH30H 


241.791 


5f,A+-4nA+ 


1.20 


1.83 


0.61 




CH30H 


241.904 


5+7E-4+9E 


0.35 


1.71 


0.19 




CH30H 


338.345 


7_iE-6_iE 


1.14 


1.43 


0.75 




CH30H 


338.404 


76E-66E 


1.85 


1.86 


0.94 




CH30H 


338.722 


7±2E-6±2E 


0.74 


2.03 


0.34 




CN 


340.248 


3-2 


23.80 


2.18 


10.3 




HNC 


362.630 


4-3 


3.83 


1.57 


2.29 




HCN 


354.505 


4-3 


41.70 


2.39 


16.4 




H"CN 


345.340 


4-3 


2.11 


1.47 


1.35 




C^^'S 


241.016 


5-4 


7.32 


1.65 


4.16 






337.396 


7-6 


5.74 


1.28 


4.21 




so 


346.528 


8,9-7,8 


4.82 


1.64 


2.77 


SK1-0MC3 


CO 


345.796 


3-2 


404.00 


- 


- 




"CO 


330.588 


3-2 


5.33 


1.12 


4.49 










69.50 


2.49 


26.2 




c"o 


337.061 


3-2 


3.31 


1.60 


1.94 




H2C0 


218.222 


3o3— 2q2 


0.57 


1.45 


0.37 




H2C0 


218.475 


322-221 


<0.20 


0.49 


0.37 




CN 


340.248 


3-2 


0.94 


3.49 


0.25 




HNC 


362.630 


4-3 


0.23 


1.25 


0.17 




HCN 


354.505 


4-3 


0.43 


1.24 


0.32 










1.74 


14.6 


0.11 
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Table 4. Physical Conditions Discerned from Carbon- 
Monoxide 





CO 




"CO 




c"o 










Source 






^peak 


TkV 


^peak 




Af(CO) 


^mod 


^obs / '^mod 




(K) 


(Kkms'') 


(K) 


(Kkms"') 


(K) 


(Kkms-') 


(10"*cm-2) 


"CO 




MMS6 


23. 


93. 


21. 


45. 


3.6 


4.8 


4.6 


2.1 


2.4 


MMS9 


28. 


225 


23. 


59. 


2.1 


3.0 


2.7 


1.3 


1.5 


FIR4 


42. 


442 


28. 


101 


2.7 


4.1 


3.4 


1.8 


1.3 


NGC2071 


67. 


770 






2.2 


7.3 


5.7 


1.5 


1.0 


PDRl 


165 


535 


105 


220 


7.2 


10. 


11. 


2.0 


1.1 


SKI 


52. 


404 


26. 


75. 


1.9 


3.3 


2.8 


1.5 


0.9 
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Table 5. Measured Total Line Intensity within the 850 /im Band 






and Contribution to Continuum Flux 








Source 


CO Line Strength!"! 


Total Line Strength™ 


Contamination 


yV(H2) 


[CO]/[H2] 




(Kknis"') 


(Kknis-'j 


% 


(10" cnr'j 


(i(r') 


MMS6 


143 


170 


1.5 


104 


0.04 


MMS9 


287 


301 


8.2 


14 


0.19 


MMS7 








22 




FIR4 


547 


934 


8.0 


47 


0.07 


NGC2071 


777 


9451^1 


8.3 


62 


0.09 


PDRl 


765 


853 


11.1 


10 


1.10 


SKI 


483 


487 


63. 


6 


0.47 



^ CO Une strengths include contribution from CO, '^CO, and C''0 

J = 3-2 transitions. 
* Total line strength includes all measured lines within the 850jum 

passband. 

The total line strength for NGC2071 does not include contribu- 
tions from methanol. Comparison with FIR4 suggests that an ad- 
ditional few percent of line contamination might be produced by 
these unobserved lines. 
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Table 6. Physical Conditions and Abundances Discerned from 
Formaldehyde 



Source 




33 


53 


T 

(K) 


n 


A'(H2CO) 
(10" 


[Hocoi/rcoi 

(10-5) 


(10-'°) 


MMS6 


n 


> 5 


0.6 


<50 


1-3 


4.0 






MMS6 


b 


> 5 


0.6 


< 50 


1-3 


4.0 






J.VJ_LV±iJ\J 












8.0 


1.7 


0.8 


MMS9 


t 


10 


0.3 


30 


1-3 


2.2 


0.8 


1.6 


MMS7 


t 


> 7 




< 40 




2.0 




0.9 


FIR4 


n 


3.3 


0.7 


90 


0.5-2 


12.0 






FIR4 




3.4 


1.8 


80 


2-10 


20.0 






FIR4 


I 










32.0 


9.4 


6.5 


NGC2071 


n 


3.2 


0.9 


90 


1-3 


15.0 






NGC2071 


b 


3.8 


0.8 


70 


1-3 


17.0 






NGC2071 


t 










32.0 


4.3 


5.1 


PDRl 


n 


1.8 




> 200 




0.8 






PDRl 


b 


2.8 


0.7 


120 


0.5-2 


6.0 






PDRl 


t 










6.8 


0.6 


6.8 



" The letter designations refer to the velocity components with nar- 
row (n), broad (b), and total (t). 

* Rsi is the ratio of the 3o3-2o2 transition versus the 322-221 transi- 
tion and R53 is the ratio of the 5o5^o4 transition versus the 3o3-2o2 
transition. 
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Table 7. Physical Conditions and Abundances Discerned from 
Methanol 



Source 


C- ' 


T 


n 


A'lCHjOH) 


[CH30H]/[C0] 


[CH30H]/[H2] 








1 1 0' cni"-' ) 


(10" rnr- ) 


( 10""^ j 




MMS6 


t 


50 


1.0 


3.0 


0.7 


0.3 


MMS9 


t 


< 30 


0.5 


1.5 


0.6 


1.1 


FIR4 


n 


70 


> 10 


30 




6.1 


FIR4 


b 


80 


5 


95 




19 


FIR4 


t 






125 


37 


25 


NGC2071 


t 


60 


0.5 


30 


5.6 


5.1 


PDRl 


t 


> 100 


5.0 


6.0 


0.6 


5.8 



The letter designations refer to the velocity components with nar- 
row (n), broad (b), and total (t). 
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Table 8. Abundances Discerned from HCN, HNC, and CN 



Source 


[HCN]/[HNC] 


Af(HCN) 


Af(CN) 


[HCN]/[CN] 


[HCN]/[H2] 


[CN(total)]/[H2] 






(lO'"* cm-2) 


(lO"* cm-2) 


(10°) 


(lO-'O) 


(lO-'O) 


MMS6 


6.3 


0.6 


<2.0 


>0.3 


0.6 


<2.4 


MMS9 


5.4 


1.3 


2.8 


0.5 


9.3 


32. 


FIR4 


11. 


8.0 


6.0 


1.3 


16. 


30. 


NGC2071 


6.0 


6.4 


19. 


0.3 


10. 


41. 


PDRl 


33. 


4.0 


2.0 


2.0 


40. 


72. 
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Table 9. Physical Conditions and Abundances Discerned from 
CS and SO 



Source 


yvccs) 


[CS]/[H2] 


Tcs 


A'(SO) 


[SO]/[H2] 


[SOl/fCS] 




(lo^ciir-) 


(K)-'") 


(K) 


(10" cnr-j 


(ur'") 


(10"j 


MMS6 


3.5 


0.3 


>50 


<2.6 


<2.5 


<7.4 


MMS9 


1.5 


1.1 










FIR4 


14. 


2.9 


>100 


<4.0 


<8.2 


<2.9 


NGC2071 


19. 


3.1 


>50 








PDRl 


29. 


29. 


>50 


0.8 


8.0 


0.3 
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5:35:00 



R.A. (2000) 



Fig. 1. The Integral Shaped Filament (ISF) at the north end 
of the Orion A molecular cloud, observed in dust continuum 
emission at 850 microns. The grey-scale ranges from to 5 
Jys with a logarithmic stretch. The positions of the six sources 
located in the ISF are indicated by arrows pointing at the peak 
continuum flux locations. Details of the sources are found in 
Tabled 
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Fig. 2. The CO J = 3-2 line profiles. Note the self-absorption 
in all sources except PDRl, and the extended line wings 
(dashed lines) in all sources except MMS6, and PDRl. 




Fig. 3. A selection of methanol lines observable with a single 
spectral setting, along with an SO Nj=S9 1% line. 
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Fig. 4. The HCN and HNC J = 4-3 transitions. 
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Fig. 5. Determination of the equilibrium physical properties of 
molecular gas from formaldehyde Une ratios. The top left panel 

presents the expected line ratio for the 3o3 - 2o2 and 322 - 22i 
transitions as functions of temperature and density (solid con- 
tours) against the expected Une ratio for the 5o5 - 4o4 versus the 
3o3 - 2q2 transitions (dashed contours). The remaining panels 
show the range of physical properties obtainable for the sources 
in the present study, with the solid lines denoting the 3o3 - 2o2 
versus 322 - 22 1 ratio and the dashed lines providing the loca- 
tion of the 5o5 - 4o4 versus 3o3 - 2o2 ratio. Two sets of dashed 
lines are shown to account for possible beam dilution (see text). 
Also, the dotted lines provide an indication of the range of un- 
certainty on the calculations. 
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CH30H 
T„t = 24.3 +/- 1.5 K 
Njit = 3.53E+14 +/- 8.17E+13 cm"^ 



(a) 



50 100 

E„p,„ / k (K) 



CH30H 
Tfit = 41.8 +/- 3.2 K 
N,„ = 3.65E+14 +/- 5-44E+13 cm"'^ 




E„pp.r / k (K) 



CH30M 
Tfit = 80-3 +/- 11.9 K 
Nfit = 2.83E+14 +/- 7.91E+13 cm"^ 




Fig. 6. Rotation diagrams for methanol. The input data are 
taken from the result of RADEX equilibrium calculations at a 
temperature of 100 K. (a) Calculation with n - 10^ cm"^. The 
slope of the data at high energies is significantly afi'ected by 
the sub-thermal excitation, and the derived temperature from 
the rotation diagram does not correspond to the physical con- 
dition. The low energy lines are much closer to being ther- 
malized, however, and the colunm density derived is close to 
the input column density of 10'"*cm"^. (b) Calculation with 
n = lO'^cm The slope of the data at high energies is still 
affected by the sub-thermal excitation, and the derived temper- 
ature from the rotation diagram continues to underestimate the 
physical condition, (c) Calculation with n = 10'^ cm"^. At this 
density the rotation diagram produces a reasonable fit to the 
input physical conditions. 
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Fig. 7. Chi-squared confidence fits to the physical parameters 
of sources using the observed methanol lines and the radiative 
transfer code RADEX. The sources are (Top Left) FIR4 narrow 
component, (Top Center) FIR4 broad component, (Top Right) 
NGC2071, (Bottom Left) MMS6, (Bottom Center) MMS9, 
(Bottom Right) PDRL Note that there is often a strong degen- 
eracy between the best fit temperature and density due to the 
lack of thermalization of the higher energy transitions. 
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Fig. 8. The range in abundances of the various molecular 
species investigated against T(H2CO). The abundances have 
been derived assuming a constant temperature and density en- 
velope. Diamonds represent sub-millimeter sources, stars rep- 
resent infrared sources, and the cross represents the PDR, 
with symbol size approximating the order unity uncertainty ex- 
pected in the measurements. For all protostellar sources the CO 
abundance is depleted by about an order of magnitude. For all 
species the infrared sources show a remarkable similarity in 
abundances, however, in the case of methanol the very broad 
> lOkms ' component from FIR4 has not been included. For 
most molecules the sub-millimeter sources have lower abun- 
dances while the PDR has the highest abundance. 
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Appendix A: Contamination of sub-millimeter 
continuum by molecular lines 

Comparison of 850//m continuum against CO 3-2 measure- 
ments in the spiral arms of external disk galaxies has revealed 
that as much as 50% of the emission within the SCUBA band- 
pass is contamination due to this single line (Tilanus, van der 
Werf, & Israel 2000). This result is not surprising since the en- 
ergy radiated in an optically thin CO j-i line by a column of 
gas with gas temperature Tg and number density n is given by 

Sji = Ay, A^(CO)Xj(Tg,n) erg s^'cm^^sr-', (A.l) 

where the column of CO is A'(CO) = x(CO) N(H2), and the 
fractional population of the jth level is Xj. For the CO 3-2 line, 
V32 = 3.458 X 10" S-' and A32 = 2.48 x 10"'' s"'. Alternatively, 
for optically thin dust the energy radiated across a passband 
centered at and with an effective width Av is 



where possible. However, as long as the line is not significantly 
broadened, the CO contamination should be minor Freeze-out 
of the CO onto dust grains decreases the importance of CO 
contamination for the densest and coldest regions. 



By.iTd) Nh, M niH Ky, Av 



erg s 'cm ^sr ', 



(A.2) 



where is the emissivity of the dust at Vd and ju is the mean 
molecular weight of the gas. Since the dust radiation in molec- 
ular clouds is usually near the Rayleigh- Jeans limit, one can 
replace the Planck function By(Td) with 2 k Tdfv{Td)/^l, where 
fv{Td) accounts for the difference between the Planck law and 
the Rayleigh- Jeans limit at frequency v. 

Thus, in the optically thin limit the ratio of the line to con- 
tinuum is independent of the column of hydrogen, depend- 
ing only on the temperature and density within the region, the 
abundance of CO in the gas phase, and the dust emissivity. 



R 



Sn\kTdl \Av 



fy,(Td)x(CO)X,(Tg,nUA.3) 



Computing the ratio of the CO 3-2 line against the JCMT 



850yum continuum band, where Av = 3 x 10"^ s yields 



7?850 - 28 



\20K/ 



x(CO) 



0.02cm2/g/ \lxlO 



h5Q(Td)X3(Tg,n).{AA) 



The ratio depends strongly on the population of the upper CO 
energy level X3 which is a function of both the gas temper- 
ature and the density. Tables lA.ll and IA.2I show X3 and R^^q 
for a range of internal gas conditions, assuming that the gas 
dust temperatures are equal. The values shown in Table IXITI 
only apply in the optically thin limit for the CO 3-2 transition 
(A^(H2) < 10^' cm"^). At higher columns, as are expected deep 
within molecular clouds, the continuum emission continues to 
increase linearly while the line emission is severely damped 
and the contamination is expected to be less influential. For ex- 
tragalactic studies, where the spatial resolution does not permit 
the separation of thick and thin portions of the cloud, it is not 
surprising that the CO emission produces a significant contam- 
ination. 

The degree to which the CO line contamination affects 
sub-millimeter continuum measurements at high optical depth 
within molecular clouds, as investigated in this paper, depends 
sensitively on the width of the CO 3-2 line and the strength of 
the extended line wings, and should be measured independently 
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Table A.l. Fractional population of the 3rd CO energy level for 
different physical conditions 





Density (cm 






Temperature (K) 


10' 


10^ 


10' 


10 


2.51e-3 


2.23e-2 


5.51e-2 


20 


1.23e-2 


9.02e-2 


1.68e-l 


30 


2.21e-2 


1.41e-l 


2.13e-l 


50 


3.66e-2 


1.95e-l 


2.24e-l 


100 


5.61e-2 


2.38e-l 


1.91e-l 


500 


8.44e-2 


2.53e-l 


1.13e-l 



Table A.2. Ratio of CO line strength to 850 yum continuum with 
SCUBA for different physical conditions, assuming [CO]/[H2] 
= 10-'' and/csso = 0.02cm2g-i. 





Density (cm ^) 






Temperature (K) 


103 




105 


10 


0.4 


3.3 


8.2 


20 


0.8 


4.0 


7.5 


30 


0.6 


3.6 


5.4 


50 


0.5 


2.6 


3.0 


100 


0.4 


1.5 


1.2 


500 


0.1 


0.3 


0.1 



